Time and distance constraints on accelerating cosmological models 
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The absence of guidance from fundamental physics about the mechanism behind cosmic acceler- 
ation has given rise to a number of alternative cosmological scenarios. These are based either on 
modifications of general relativistic gravitation theory on large scales or on the existence of new fields 
in Nature. In this paper we investigate the observational viability of some accelerating cosmological 
models in light of 32 age measurements of passively evolving galaxies as a function of redshift and 
recent estimates of the product of the cosmic microwave background acoustic scale and the baryonic 
acoustic oscillation peak scale. By using information-criteria model selection, we select the best-fit 
models and rank the alternative scenarios. We show that some of these models may provide a better 
fit to the data than does the current standard cosmological constant dominated (ACDM) model. 



PACS numbers: 98.80.-k, 95.36. 



I. INTRODUCTION 



98.80.Es 



There is fairly convincing observational evidence that 
the Universe is now undergoing accelerated expansion. In 
the context of Einstein's general theory of relativity (GR) 
this amounts to saying that some sort of dark energy, 
constant or that varies slowly with time and in space, 
dominates the current energy budget of the cosmos. On 
the other hand, it is also conceivable that the current 
accelerated expansion could be an indication that GR 
does not provide an adequate description of gravity on 
cosmological length scales. For recent reviews, see 

While there are many options, if GR is assumed to be 
the correct model of gravitation, the current consensus 
opinion is that the simplest possible driver of this ac- 
celerated expansion is a cosmological constant, A, a spa- 
tially homogeneous fluid with equation of state p\ = —p\ 
(where p\ and p\ are the fluid pressure and energy den- 
sity) or equation of state parameter u> = — 1. This "stan- 
dard" ACDM model @, where cold dark matter (CDM) 
is the second largest contributor to the cosmological en- 
ergy budget after A, allows for a good fit to many curent 
observations. 1 However, the needed A energy density 
scale of order 10 -3 eV is very small and difficult to rec- 
oncile with the much higher value that a naive appli- 
cation of quantum mechanics suggests. Furthermore, as 
the density of nonrelativistic matter (the CDM and bary- 
onic matter) decreases during the cosmological expansion 
while the A energy density remains constant, these two 
energy densities are comparable in magnitude only for a 
relatively short time and it is puzzling why we happen to 
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1 The "standard" CDM model of structure formation might be 
inconsistent with some observations (see, e.g., Q). 



observe near this epoch (this is the so-called coincidence 
puzzle). 

Motivated by the fact that instead of a cosmological 
constant, if the energy density responsible for the accel- 
erated expansion — the dark energy — also decreased in 
time, the energy densities of nonrelativistic matter and 
dark energy would be comparable for a somewhat longer 
period, it is of interest to consider models where the dark 
energy density decreases with time. A simple example is 
the XCDM parametrization of dark energy, where dark 
energy is modeled by a spatially homogeneous fluid with 
equation of state px = ^xPx with tax < 0. While the 
XCDM parametrization is widely used, it is well known 
that it is not a complete parametrization as it cannot 
describe spatial inhomogeneities (see, e.g. Q). 

The ifiCDM model 5], in which dark energy is taken 
to be a scalar field <j) with a potential energy density 
V{<j>) that results in a dark energy density that is close 
to spatially homogeneous but slowly decreasing in time, 
provides a complete description of dark energy, unlike 
the XCDM parametrization. See, e.g. [6| for more recent 
discussions about this model. 

Rather than introducing a new dark energy compo- 
nent, explanation of the accelerated cosmological expan- 
sion has been sought in modifications of GR, including 
the possible existence of extra dimensions (e.g. [2|) or by 
modifying GR by adding terms proportional to powers 
of the Ricci scalar R to the Einstein-Hilbert Lagrangian 
(e.g. 1). 

Given the state of uncertainty that remains over 
whether A may or may not provide a complete expla- 
nation for the observed cosmic acceleration, a very im- 
portant way to improve our understanding of this phe- 
nomenon is to use cosmological observations to constrain 
mechanisms of cosmic acceleration. In this regard, it 
is worth emphasizing that using different observational 
techniques, as well as combinations of them, is particu- 
larly important for a more reliable determination of cos- 
mological parameters, since different methods may rule 
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out different regions of the parameter space and, there- 
fore, may be complementary to each other. 

In this paper we combine distance data with low and 
high-redshift time measurements to constrain acceler- 
ating cosmologies. Specifically, we use age measure- 
ments of 32 passively evolving galaxies Q (in the range 
0.117 < z < 1.845) to perform a statistical analysis of 
the age-z test. In order to better constrain the parameter 
spaces of these models, we combine the age- 2 data with 
a recent estimate of the ratio of the CMB acoustic scale 
£a and the baryonic acoustic oscillation (BAO) peak, the 
so-called CMB/BAO ratio [l(|. We show that for some 
models the age-z plus CMB /BAO ratio data analysis re- 
sults in restrictive joint constraints on these cosmologies. 
We discuss our results in terms of two information cri- 
teria, namely, the Akaike Information Criterion [ll"j | and 
the Bayesian Information Criterion [l2[ . 

The outline of our paper is as follows. In Sec. 2 we 
present the basic equations of, and main motivations for, 
the accelerating cosmological models we consider in our 
analyses. The age-redshift test and the joint analysis of 
this data with the CMB/BAO ratio data are outlined in 
Sec. 3. Our results are presented and discussed in Sec. 4. 
In Sec. 5 we summarize our main results. 



II. NONSTANDARD ACCELERATING 
COSMOLOGIES 

As mentioned in the previous Section, theoretical mod- 
els of cosmic acceleration can be broadly classed into two 
groups: those that introduce a new component to the 
composition of the Universe (this includes the "standard" 
ACDM case) and those that invoke a modification to the 
equations governing gravity. In what follows, we briefly 
discuss some of the more popular nonstandard models 2 
and examine whether they are consistent with current 
age and distance data. Unless specified otherwise, we 
assume vanishing spatial curvature. 



A. Quintessence scalar field 

For some potential energy density functions, V(4>), the 
</>CDM model partially resolves two theoretical puzzles of 
the ACDM model. For example, with an inverse power 
potential V(<p) oc <j)~ a with a > [5] the solution of the 
coupled nonlinear equations of motion is an attractor for 
which the scalar field energy density, p^, decreases less 
rapidly than the dominant component of the energy den- 
sity in the radiation- and matter-domianted epochs and 
so eventually comes to dominate, thus partially resolv- 
ing the coincidence puzzle. Also, the slow decrease in 
time of the dark energy density in this model results in 



a very small dark energy density now because the Uni- 
verse is old. Consequently, instead of a new fundamental 
energy density scale of order a meV as in the ACDM 
model, the 0CDM model requires a higher fundamen- 
tal energy density scale that is determined by when the 
phenomenological 0CDM model first provides an appro- 
priate description of cosmology. 

With V(<fi) — n(j)~ a , where k is a constant, the dy- 
namics of the tfiCDM model is obtained by solving the 
equations 



E l = 



8irG 



(p m + P<p) , 



(1) 



(2) 



where H is the Hubble parameter, p m is the pressureless 
matter (baryons and cold dark matter) energy density 
and a dot denotes a derivative with respect to time. The 
parameter a describes the steepness of the scalar field 
potential, with larger (smaller) values of a corresponding 
to faster (slower) evolution of the scalar field. In the limit 
a — >• the spatially-flat </>CDM model reduces to the 
spatially-flat ACDM case. The observational viability 
of this model has been extensively investigated in the 
literature and we refer the reader to 11411 . 



B. Chaplygin gas 

From the cosmological viewpoint, the main distinction 
between pressureless CDM and dark energy is that the 
former agglomerates at small length scales whereas the 
latter is a smooth component on these scales. Recently, 
the idea of a unified description for CDM and dark energy 
has received much attention (see, e.g., |l5()- An interest- 
ing attempt in this direction was suggested in Ref. [l6[ 
and further developed in Ref. [TtJ • It uses to an exotic 
fluid, the so-called Chaplygin gas (Cg), whose equation 
of state is given by 



PCg 



A 



(3) 



Actually, this equation for a / 1 is a generalization of 
the original Cg equation of state and was originally pro- 
posed in [l8| (see also [Ii|). In our analyses, we will con- 
sider both the original Cg (a — 1) and its generalization 
(GCg). _ 

Inserting Eq. §3§ into the energy conservation equation 
gives the expression for the GCg energy density 



PCg = PCg 



A s + (1-A S 



-3(H 



(4) 



For other alternative models, see US 



where a(t) is the cosmological scale factor and A s — 
A/ Peg" 1S a quantity related to the sound speed for the 
generalized Chaplygin gas today (throughout this pa- 
per the subscript "0" denotes present-day quantities). 
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As can be seen from the above equations, the Chap- 
lygin gas interpolates between epochs dominated by 
non-relativistic matter [pcg{a — >• 0) ~ const. /a 3 ] and 
by a negative-pressure time- independent dark "energy" 
[pc g (a -> oo) ~ y/A\. 



C. Coupled quintessence 

A phenomenological attempt at alleviating the cosmo- 
logical constant problems mentioned earlier is based on 
allowing the dark matter and dark energy to interact. 
For this class of models, the energy balance equation be- 
tween a quintessence component <p an d CDM particles 
must be of the form 



Pcdm ~t~ 3 Pcdm 

a 



3-/9^(1 + w$) . (5) 



The dark matter-dark energy interaction implies that the 
energy density of the former component must dilute at a 
different rate compared to its standard evolution p c dm oc 
a~ 3 (Alcaniz 2006). If we characterize this non-standard 
evolution by an arbitrary function e(a) such that p c dm = 
pcdm,oa>~ 3+e ^ the energy density of the <fi component 
can be written as (Costa 2010) 



P<t> 



P4>,0 + Pcdm 



f 1 Ha 



e(o) + ae Ind] 



da 



-3(l+u^) 



-3ro^-e(a) 

(6) 

For Wcf, = —1 the above expression reduces to the 
vacuum decay model discussed by Costa & Alcaniz 
(2010), whereas for a constant e one recovers the cou- 
pled quintessence scenario studied by Jesus ct al. (2008). 
Note also that the (incomplete) XCDM parametrization 
is recovered for e = 0. 

For simplicity, we restrict our analyses to coupled 
quintessence models in which e = eo = const, and 
Wff, — —1, so that 



eoPcdmfi -3+e 

3 - e 



(7) 



which is mathematically equivalent to the decaying vac- 
uum scenarios proposed by Wang & Meng (2005) and 
Alcaniz & Lima (2005) (see Amendola 2000; Dalai et al. 
2001; Wu et al. 2008; Carneiro et al. 2008, for other sce- 
narios of coupled quintessence). Since we are assuming 
interaction only between the components of the dark sec- 
tor, in our analyses of this model we fix the baryonic mat- 
ter fractional energy density contribution at = 0.0416, 
as given by current CMB anistropy measurements (Ko- 
matsu et al. 2009) for Hubble constant H = 74.2 km s _1 
Mpc" 1 . 



mcnt to the energy content of the Universe. 3 This mo- 
tivates studies of cosmological models with extra spatial 
dimensions that might be able to explain the accelera- 
tion as well as provide a possible rationale for the huge 
difference between the electroweak and Planck scales 
mpi/niEw ~ 10 16 (Arkani-Hamed et al. 1998; Randall 
& Sundrum 1999). 

Many extra-dimensional models have been discussed. 
Here we consider the DGP model (Dvali et al. 2000), 
a self-accelerating 5-dimensional model with a non- 
compact, infinite- volume extra dimension, whose grav- 
itational dynamics is governed by a competition between 
a 4-dimensional Ricci scalar term induced on the brane 
and an ordinary 5-dimensional Einstein-Hilbert action 
(for a review of DGP phenomenology, see Lue 2006). 
The Friedmann equation derived from this gravitational 
action take the form 



H' + — = 

a* 



2,M 2 pl + 4r2 



1 

2^ 



(8) 



where p is the energy density of the cosmic fluid, k is the 
spatial curvature, and r c = Mj> ; /2M| is the crossover 
scale defining the gravitational interaction among parti- 
cles located on the brane. From the above equation we 
can also define the density parameter associated with the 
crossover radius, fl rc — \/Ar 2 H^ (see also Dcffayct et al. 
2001; Alcaniz 2002; Guo et al. 2006). For this scenario 
we analyze both flat (flat DGP) and arbitrary-curvature 
cases (DGP). 



E. Cardassian expansion 

Also based on extra dimensions, Freese & Lewis (2002) 
proposed the Cardassian model in which the Universe is 
spatially flat, matter dominated, and currently undergo- 
ing accelerated expansion. Such a behavior is obtained 
from a modified Friedmann equation on our observable 
brane with H 2 — g(p m ), where g(p m ) is an arbitrary 
function of the matter energy density. Although com- 
pletely different physically, the original Cardassian model 
with g(pm) = A m + B n (where A and B are functions of 
the matter density p m and n and m are free parameters 
of the model) predicts the same observational effects as 
the XCDM parametrization for cosmological tests based 
on the evolution of the Hubble parameter with redshift. 

Here, we consider the generalized Cardassian (GCard) 
model described by Wang et al. (2003) whose Friedmann 



D. DGP 

A different possible explanation of cosmic accelera- 
tion is a modification in gravity instead of an adjust- 



3 As mentioned earlier, in addition to the extra dimensional models 
we consider, another way of modifying gravity is by adding terms 
proportional to powers of the Ricci scalar R to the Einstein- 
Hilbert Lagrangian, the so-called f(R) gravity model. For some 
recent applications of f(R) gravity in cosmology see Sotiriou 
(2009), de Felice & Tsujikawa (2010), and references therein. 
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equation is given by 



H- 



&irGp ri 



\ Pm J 



9(1-™)' 



1/'/ 



(9) 



where n and q are two free parameters to be adjusted to 
fit the data and p C ard = p m ,o( l + z Card) 3 is the energy 
density at which the two terms inside the bracket become 
equal. Note that, although matter dominated, GCard 
models are accelerating and may still reconcile current 
indications for a spatially-flat Universe from CMB obser- 
vations with clustering estimates that point consistently 
to fl m ~ 0.3. In these scenarios, this discrepancy is ex- 
plained through a redefinition of the value of the critical 
density from Eq. © (see, e.g. Freese & Lewis 2002; Al- 
caniz et al. 2005 for discussions). 
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III. COSMOLOGICAL CONSTRAINTS 

A. The age-redshift test 

The theoretical age-z relation [t(zi)] of an object at 
redshift Zi can be written as (Sandage 1988) 



t( Zl ,p) = / 
Jo 



dx' 



x'H(x',p) ' 



(10) 



where p stands for the parameters of the cosmological 
model under consideration and H(x, p) is the normalized 
Hubble parameter. From the observational viewpoint, 
the total age of a given object (e.g., galaxies) at redshift 
z is given by £ obs (^) = ta(zi) + r, where where tc(zi) 
is the estimated age of its oldest stellar population and 
r is the incubation time or delay factor, which accounts 
for our ignorance about the amount of time since the 
beginning of structure formation in the Universe until 
the formation time of the object of interest. 

To perform our age-z analyses in the next section, we 
use age estimates of 32 old passive galaxies distributed 
over the redshift interval 0.117 < z < 1.845 (Simon et 
al. 2005) as listed in Table 1 of Samushia et al. (2010), 
and assume a 12% one standard deviation uncertainity 
on the age measurements (R. Jimenez, private commu- 
nication 2007). The total sample is composed of three 
sub-samples: 10 field early-type galaxies from Treu et al. 
(1999; 2001; 2002), whose ages were obtained by using 
the SPEED models of Jimenez et al. (2004); 20 red galax- 
ies from the publicly released Gemini Deep Deep Survey 
(GDDS), whose integrated light is fully dominated by 
evolved stars (Abraham et al. 2004; McCarthy et al. 
2004); and the 2 radio galaxies LBDS 53W091 and LBDS 
53W069 (Dunlop et al. 1996; Spinrad et al. 1997; Nolan 
et al. 2001). As discussed by McCarthy et al. (2004), 
the GDDS data seem to indicate that the most likely star 
formation history is that of a single burst of duration less 
than 0.1 Gyr, although in some cases the duration of the 
burst is consistent with Gyr, which means that the 



FIG. 1: Contours of A X l = 6.17 (2cr) for the ACDM model 
(left panel) and the spatially-flat XCDM parametrization 
(right panel). The contours correspond to the joint analy- 
sis involving age-z plus CMB/BAO data. In the left panel 
(ACDM) the thin orange diagonal line demarcates spatially- 
flat models while in the left panel (XCDM) the thin orange 
horizontal line indicates models with a time-independent cos- 
mological constant. 



galaxies have been evolving passively since their initial 
burst of star formation. 

The likelihood function C oc [exp — Xa ffe (z; p, r)/2] is 
determined by a \ 2 statistics 



<l<- 



(H 0: r,p)=J2 



[t(zi,p) - t G (zi) - t\ 



(11) 



where of . stands for the uncertainties on the age mea- 
surements of each galaxy in our sample. An important 
aspect in our analysis that is worth emphasizing con- 
cerns the delay factor r. Note that in principle there 
must be a different Ti for each object in the sample since 
galaxies form at different epochs. Here, however, since 
we do not know the formation redshift for each object, 
we assume a uniform delay factor r that we treat as a 
"nuisance" parameter and marginalize over it. Similarly, 
we also marginalize over the present value of the Hubble 
parameter Hq. 4 Note also that, in order to avoid dou- 
ble counting of information with CMB data discussed in 
the next section, we have not included any prior on the 



4 A variant of this test uses both measurements of the Hubble 
parameter as a function of redshift (see, e.g. Samushia & Ratra 
2006; Samushia et al. 2007; Dev et al. 2008; Gong et al. 2008; 
Fernandez-Martinez & Verde 2008; Wei 2009; Yang & Zhang 
2009) and lookback time measurements built from estimates of 
the total age of the Universe (Capozziello et al. (2004), Pires et 
al. (2006), Dantas et al. (2007, 2009), Samushia et al. (2010) 
and Samushia (2009). 
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FIG. 2: Contours of \t corresponding to the joint analysis 
discussed in the text for two of the accelerating cosmologies. 
Left panel shows the spatially-flat 0CDM model with V{4>) oc 
4>~ a while the right panel is for the Generalized Chaplygin 
gas case. The contours are drawn for A% 2 = 6.17 (2cr). For 
a — the spatially-flat (pCDM model reduces to the spatially- 
flat ACDM case. 

age of the Universe from CMB data, as usually done in 
lookback time tests. 



B. CMB/BAO ratio 

The two major inputs involving acoustic oscillations 
come from the CMB and baryon oscillations data. Al- 
though the CMB shift parameter TZ and the BAO param- 
eter A have been commonly used to constrain nonstan- 
dard models, their use may not always be entirely appro- 
priate because the values used were obtained in the con- 
text of an extended XCDM parametrization, which can 
be considered a good approximation only for some classes 
of dark energy models (see, e.g. Doran et al. 2007). 5 

Here, we follow Sollerman et al. (2009) and use a more 
model-independent constraint derived from the product 
of the CMB acoustic scale I a = Trd,A{z st )/r s {z^) and the 
measurement of the ratio of the sound horizon scale at the 
drag epoch and the BAO dilation scale, r s (zd)/Dy (zbao) 
[<iA(z*) is the comoving angular-diameter distance to 
recombination and r s (z*) is the comoving sound hori- 
zon at decoupling]. By combining the ratio r s (zd — 
1020)/r s (z* = 1090) = 1.044 ± 0.019 (Komatsu et al. 



5 For the BAO parameter A, for instance, it is implicitly assumed 
that the evolution of matter density perturbations during the 
matter-dominated era must be similar to the ACDM case and 
also that the comoving distance to the horizon at the time of 
equilibrium between matter and radiation energy densities must 
scale with (S7 m iig) -1 . 



2009) with the measurements of r s (zd)/ Dv(zbao) a t 
zbao = 0.20 and 0.35 from Percival et al. (2009), Soller- 
man et al. (2009) found (with one standard deviation 
error bars) 

d A (z*)/D v (0.2) = 17.55 ± 0.65 

d A {z*)/D v {0.35) = 10.10 ±0.38 , 

which we use in our analyses along with the LT-z data 
discussed earlier. We account for the correlations in the 
measurements by following Sollerman et al. (2009). 

IV. RESULTS 

Figures 1, 2, and 3 show some of the main results of 
our analyses. In these figures, we show 2a contours for 
the the joint analysis involving age-z plus CMB/BAO 
data. Note that, although not providing very restric- 
tive constraints in all the cases analized, these time and 
distance cosmological data can place tighter bounds on 
some cosmological parameter spaces. In particular, for 
the ACDM (Fig. 1 left panel) and DGP (Fig. 3 center 
panel) cases the x 2 contours for age-z and CMB/BAO 
data are roughly orthogonal, which results in tight cos- 
mological constraints when used together (see, e.g., Pires 
et al. 2006; Samushia et al. 2010). 

The shadowed area in the left panel of Fig. 3 for the 
coupled quintessence case corresponds to the thermody- 
namical constraint on the interaction parameter (eo > 0) 
discussed by Alcaniz & Lima (2005). In agreement with 
other recent analyses (Bertolami et al. 2007; Costa et 
al. 2009), we note that half of the 2a interval for eo lies 
in the negative region of the plot, which means that an 
energy flow from dark energy to dark matter is allowed 
from this combination of data (see Eq. ([7])). 

In order to select the best-fit models, we use two infor- 
mation criteria, namely, the Bayesian information crite- 
rion (BIC) and the Akaike information criterion (AIC): 

BIC = -21n£ + fclnA , (12) 

AIC = -21n£ + 2k , (13) 

where C is the maximum likelihood, k is the number 
of model parameters, and TV is the number of data 
points used in the fit. We refer the reader to Lid- 
die (2004), Godlowski & Szydylowski (2005), Biesiada 
(2007), and Davis et al. (2007) for reviews on the back- 
ground for the use of these ICs. ICs provide an interest- 
ing way to obtain a relative ranking of the observational 
viability of different candidate models. Thus, the im- 
portant quantities in these analyses are the differences 
AAICi = AIC, - AIC min and ABIC; = BIC, - BIC min 
calculated over the whole set of scenarios i = 1 , . . . , N 
with the best-fit model being the one that minimizes 
the AIC and BIC factors. Here, we follow Liddle (2004) 
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FIG. 3: Contours of constant Axt = 6.17 for the spatially-flat coupled quintessence (left panel), DGP (middle panel), and 
generalized Cardassian (right panel) models. The shadowed area in the left panel (coupled quintessence) corresponds to the 
thermodynamical constraint on the interaction parameter (eo > 0) discussed by Alcaniz & Lima (2005). The thin orange line 
in the middle panel (DGP) demarcates spatially- flat models. 



TABLE I: Summary of the age-z IC Results 



Model 



k AAIC ABIC 



TABLE II: Summary of the age-z and CMB/BAO IC Results 



Model 



fc AAIC 



ABIC 



Flat XCDM 
Flat ACDM 
ACDM 
Flat DGP 
DGP 
Flat 0CDM 
Chaplygin 
Gen. Chaplygin 
Gen. Cardasssian 
Coupled Quint. 



2 (5) 2.78 

1 (3) 0.94 

2 (8) 2.85 

1 (1) 0.00 

2 (2) 0.93 
2 (7) 2.80 

1 (4) 1.91 

2 (6) 2.79 

3 (9) 4.58 
2 (6) 2.79 



3.97 (4) 
0.94 (2) 
4.03 (7) 
0.00 (1) 
4.11 (8) 
3.99 (6) 
1.90 (3) 

3.98 (5) 
7.40 (9) 
3.98 (5) 



Flat XCDM 
Flat ACDM 
ACDM 
Flat DGP 
DGP 
Flat 0CDM 
Chaplygin 
Gen. Chaplygin 
Gen. Cardasssian 
Coupled Quint. 



2 (9) 1.87 

1 (1) 0.00 

2 (7) 1.58 

1 (5) 1.45 

2 (4) 1.40 
2 (6) 1.52 

1 (3) 1.33 

2 (8) 1.80 



3.14 (9) 
0.00 (1) 
2.85 (6) 
3.00 (7) 
2.67 (3) 
2.80 (5) 
1.54 (2) 
3.07 (8) 



3 (10) 4.23 6.54 (10) 
2 (2) 1.04 2.69 (4) 



and adopt the criterion that a difference of 2 is consid- 
ered positive evidence against the model with the higher 
BIC, while a ABIC of 6 is considered strong evidence 
against the model (we also adopt similar rules for the 
AIC). Note also that, since in our analysis N/k < 40, 
we use the version of AIC corrected for small samples 
AIC C = AIC + 2k(k + l)/(N-k - 1) (Sugiura 1978; 
Burnham & Anderson 2002, 2004; Liddle 2004). 

Table 1 shows a summary of the IC results for the age-z 
test. The best-fit model from both ICs is the spatially- 
flat DGP model, in contrast to earlier analyses, except 
for the result shown in Table 2 of Sollerman et al. (2009) 
based on the SDSS SNe la sample (MLCS-fit) and the 



CMB/BAO ratio data discussed above. Following the 
rule for the BIC discussed earlier, we note that there is 
mostly only positive evidence in favor of the best-fit sce- 
nario relative to most of the models considered (ABIC 
~ 2). The exceptions are the DGP model with arbitrary 
spatial curvature (ABIC = 4.11) and the generalized Car- 
dassian expansion scenario (ABIC = 7.40), where the ev- 
idence in favor the best-fit model is considerably stronger. 

These results arc considerably modified when the 
CMB/BAO data are added to the analysis (Table 2). 
The best-fit model now is the fiat ACDM scenario, and 
there is now strong evidence against only the generalized 
Cardassian scenario (ABIC = 6.54). It is also worth men- 
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tioning that five (Coupled Quint., DGP, Chaplygin gas, 
spatially-flat </>CDM, generalized Chaplygin gas) out of 
the seven non-standard cosmologies considered in this pa- 
per cannot be ruled out by this combination of data. The 
results found in our joint analyses reflect the complemen- 
tarity between age-z and CMB/BAO measurements and 
confirm our initial arguments for a joint analysis involv- 
ing time and distance-based cosmological observables. 

V. CONCLUSION 

We have investigated the observational viability of 
some non-standard accelerating cosmologies in light of 
current observational data. In contrast to most recent 
analyses, which are essentially based on distance mea- 
surements to a particular class of objects or physical 
rulers, we have used complementary measurements in- 
volving time (age-z) and distance (CMB/BAO) data. 

We have shown that some alternative models not only 
survive the investigation but in some cases also provide 
a better fit to the data than does the current standard 



ACDM cosmological scenario (see Tables 1 and 2). These 
results are in good agreement with recent analyses using 
SNe la plus CMB/BAO data (see, e.g., Sollerman et al. 
2009) and reinforces the need for investigations of alter- 
native mechanisms of cosmic acceleration that rely on 
different kinds of observations. 

While current age-z data alone are not able to strongly 
discriminate between different cosmological models, we 
expect a new higher quality set of high-z age measure- 
ments to soon become available (R. Jimenez, private 
communication, 2009). We anticipate that this new data 
will prove very useful in narrowing down the range of 
viable cosmological models. 
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